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ABSTRACT 



Context. In nearby clusters early-type galaxies follow isotropic orbits. The orbits of late-type galaxies are instead characterized by 
slightly radial anisotropy. Little is known about the orbits of the different populations of cluster galaxies at redshift z > 0.3. 
Aims. We investigate the redshift evolution of the orbits of cluster galaxies. 

Methods. We use two samples of galaxy clusters spanning similar (evolutionary corrected) mass ranges at different redshifts. The 
sample of low-redshift (z ~ 0.0 - 0.1) clusters is extracted from the ESO Nearby Abell Cluster Survey (ENACS) catalog. The sample 
of high-redshift (z ~ 0.4 - 0.8) clusters is mostly made of clusters from the ESO Distant Cluster Survey (EDisCS). For each of 
these samples, we solve the Jeans equation for hydrostatic equilibrium separately for two cluster galaxy populations, characterized 
by the presence and, respectively, absence of emission-lines in their spectra ('ELGs' and 'nELGs' hereafter). Using two tracers of the 
gravitational potential allows to partially break the mass-anisotropy degeneracy which plagues these kinds of analyses. 
Results. We confirm earlier results for the nearby cluster sample. The mass profile is well fitted by a Navarro, Frenk & White (NFW) 
profile with concentration c = 4. The mass profile of the distant cluster sample is also well fitted by a NFW profile, but with a 
slightly lower concentration, as predicted by cosmological simulations of cluster-sized halos. While the mass density profile becomes 
less concentrated with redshift, the number density profile of nELGs becomes more concentrated with redshift. In nearby clusters, 
the velocity anisotropy profile of nELGs is close to isotropic, while that of ELGs is increasingly radial with clustercentric radius. 
In distant clusters the projected phase-space distributions of both nELGs and ELGs are best-fitted by models with radial velocity 
anisotropy. 

Conclusions. No significant evolution is detected for the orbits of ELGs, while the orbits of nELGs evolve from radial to isotropic 
with time. We speculate that this evolution may be driven by the secular mass growth of galaxy clusters during their fast accretion 
phase. Cluster mass density profiles and their evolution with redshift are consistent with predictions for cluster-sized halos in A Cold 
Dark Matter cosmological simulations. The evolution of the nELG number density profile is opposite to that of the mass density 
profile, becoming less concentrated with time, probably a result of the transformation of ELGs into nELGs. 

Key words. Galaxies: clusters: general ~ Galaxies: kinematics and dynamics 



1. Introduction 

' The distributions of cluster early-type and late-type galaxies 
. (ETGs and LTGs hereafter) have long been known to be dif- 
' ferent ([Dressier 1980; see also Biviano 2000 for a review). Most 
striking, and hence the first to have been discovered, is the dif- 
. ference in the spatial distributions, ETGs living in higher den- 
' sity regions than LTGs, the so-called "morphology-density rela- 
' tion" (MDR hereafter). In relaxed clusters density is an almost 
monotonic decreasing function of clustercentric radius, hence 
the r elation is often describ ed as a morphology-radius relation 
(e.g. IWhitmore et all [19931) . The MDR h as been found to ex- 
ist in rich cluster s up to redshift z ~ 1 dPostman et al J [20051; 
ISmithet al.ll2005h . but the global fractions of ETGs and LTGs 
evolve with z. The fraction of ETGs in clusters decreases with 
z quite rapidly up to z ~ .5 (Dressi er et alJll997[ ; iFasano et alj 
120001; iPostman et al.l200l [Smith et al.l2005l) . then seems to flat- 
ten out to z ~ 1 dDesai et al.ll2007l) . 

Another aspect of the segregation of different galaxy pop- 
ulations in clusters is the difference in the velocity distribu- 
tions of ETGs and LTGs. In clusters, the velocity distribution 
of LTGs is broader than that of ETGs, i.e. LTGs have a larger 
line-of-sight (los hereafter) velocity dispersion (crios hereafter) 



flian E TGs (lTammannll97H iMoss & Dickenslll977t ISodre et all 
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119891 ; iBiviano et al. 1992). The difference is not only in the 
global (Xios, but also in the cr/o^-profiles, ' hotter' and steeper (at 
least in the inner regions) for the LTGs dCarlberg et alJll997a: 
IBiviano et all 1 9971; (Adami et al.lll998al) . 

While the z-evolution of the MDR has been extensively 
studied and described, less is known about the z-evolution of 
the morphological segregation in velocity space, because of the 
difficulty of obtaining large spectro scopic s a mples of cluster 
galaxies at high-z. As remarked by iBivianol (|2006|) . the crios- 
profiles of both ETGs and LTGs in z ~ 0.07 c lusters (from 
ENACS, the ESO Nearby Abell Cluster Survey, 'Katge rt et all 
il996; Biviano et al. 19 97) are remarkably similar to the profiles 
of red, and respectively, blue galaxies in z ~ 0.3 clusters (from 
CNOC, the Canadian Network for Observational Cosmology 
survev. [Carlberg et al.| [r997b). but little is known about the evo- 
lution at still higher z. 

Instead of looking separately at the spatial and velocity seg- 
regation of cluster galaxies, it is possible to use the joint infor- 
mation coming from their 2-dimensional projected phase-space 
distribution in the los velocity vs. clustercentric radius diagram. 
The projected phase-space distribution of a given class of clus- 
ter galaxies is the observable that enters th e Jeans equation for 
the equiUbriurn of a g alaxy system (see, e.g. lBinnev & Tremaind 
119871: IBivianol [2OO8I) . hence separating dififerent cluster galaxy 
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populations on the basis of their projected phase-space distri- 
butions is a way of identifying diff erent, independent t racers 
of the cluster gravitational potential. iBiviano et alj (|2002|) have 
shown that cluster ETGs and LTGs have significa n tly di fferent 
projected phase-space distributions; ^ Katgert et alJ (12004 ) have 
then used cluster ETGs to determine the average mass pro- 
file of massive clusters. In order to do that, they first had to 
constrain the velocity anisotropy profile of ETGs. They did so 
by comparing the velocit y distribution of ETGs wi th distribu- 
tion function models from lvan der Marel et al.l (12000) , and found 
that ETGs move on nearly isotropic orbits. On the other hand, 
LTGs were found to move on slightly radially anisotropic orbits, 
with an increa sing radial anisotropy at larger clustercentric radii 
(iBiviano & K atgert 2004). 

At higher redshift, in the CNOC cluster s, red galaxies were 
found to move along nearly isotropic orbits dvan der Marel et all 
I2OOOI) . similarly to ETGs in nearby ENACS clusters, and blue 
galaxies were found to be in equilibrium within the cluster po- 
tential, despite having a diff'erent projected phase-spa ce distribu- 
tion from that of red galaxies ( Carlberg et alJ[T997bl) . Although 
a solution for the velocity anisotropy of the blue CNOC cluster 
galaxies has not been derived, the similarity of their projected 
phase-space distribution to that of LTGs in E NACS suggests 
that th ey similarly move on slightly radial orbits dBivian ol l2006i 
I2OO8I) . 

Hence, while there is significant evolution in the relative 
fractions of early-type (red) and late-type (blue) cluster galax- 
ies, their orbital anisotropics do not seem to evolve over the same 
0-0.3 redshift range. Benatovet al. (2006) do claim significant 
orbital evolution for the whole cluster galaxy population, on the 
basis of the analyses of three low-z and two z ~ 0.2 - 0.3 clus- 
ters. Specifically, they find galaxies in the higher-z clusters to 
have more radially anisotropic orbits than galaxies in the lower- 
z clusters. Since they consider all cluster galaxies together in 
their analysis, and since LTGs are known to be characterized by 
radial orbital anisotropv. lBenatov et alj (l2006l) 's result could be 
explained by the increasing fraction of LTGs with z, without the 
need for any evolution of the orbits of either ETGs or LTGs. 

The lack of a significant evolution in the orbital anisotropy 
of, separately, early-type (red) and late-type (blue) cluster galax- 
ies from z ~ 0.05 to z ~ 0.3, coupled with the significant evo- 
lution in the relative fractions of these two populations over the 
same redshift range is intriguing. It suggests that the change of 
class from blue, LTG, to red, ETG, goes together with the orbital 
change, from moderately radial to isotropic. The mechanisms 
by which this evolution occurs are not known. To gain further 
insight into this evolutionary process, in this paper we extend 
the analysis of galaxy orbits in clusters to higher redshifts than 
examined so far We base our analysis on the sam ple of clus- 
ters f rom the ESO Distant Cluster Survey (EDisCS. IWhite et al.1 
l2005h . that span the redshift range ^ 0.4-LO. 

We adopt //q = 70 km s ' Mpc ', = 0.3, Qa = 0.7 
throughout this paper 

2. Two data-sets 

The high-z (distant) cluster galaxies dat a-set used in this 
paper has beeii gath ered in the EDisCS dWhite et al.l 120051: 
iPoggianti et al.ll2006l) . a survey of 20 fields containing galaxy 
clusters in the z-range 0.4-LO. Multi-band optical and near- 
infrared photome try for these fiel ds has been obtained us- 
ing VLT/F0RS2 d White etal.1 120051) . and NTT/SOFI (Aragon- 
Salamanca et al., in prep.). Imaging with the ESO WFI has 
also been obtained as well as HST/ACS mosaic images for 10 



cluster fields dPesai et al.ll2007h and deep optical spectroscopy 
for 18 cluster fields using VLT/FORS2 dHallidav et al.ll2004i: 
iMilvang- Jensen et al ]|2008l) . 

Three of the 18 ED isCS clusters (CI 110 3, CI 1119, and 
CI 1420, see Table 1 in Poggia nti et al.l 12006*) have very low- 
masses with velocity dispersion crigs< 250 km s"', more typical 
of groups than clusters. We exclude them from our data-set since 
including them would break the expected homology of cluster 
mass profiles that is a requirement for the stacking procedure 
(see Sect. O. We are thus left with 15 EDisCS clusters. In order 
to increa se our data-set we add four clusters from t he MORPHS 
data-set dDressler et al.lll999tlPoggianti et al.l 19991) with masses 
in the same range covered by the 15 EDisCS clusters. These 
cl usters (CI 0024 CI 03 03, CI 0939, and CI 1601; see Table 2 
in lPoggianti et al.l l2006') have sufficiently wide spatial coverage 
(> 0.5 r2oofl needed for the determination of the mass profile, 
and homogeneous photometry, needed for the determination of 
the radial incompleteness (see Sect.[3j. 

As a reference low-z (nearby) sample of clusters, we use th e 
59 ENACS clusters studie d in detail by IBiviano et all d2002l) ; 
IBiviano & KatgertI (|2004 ; iKatgert e^a l' (2004). A full descrip- 
tion of the ENACS can be found in [Katgert et al. ( 1996, 1998). 

We identify and reject interlopers in both the high- 
and the low-z clus ters following the procedure described in 
IBiviano et al. I (120061) . which is based on the identification of sig- 
nificant gaps in redshift space dGirardi et al.lll993l) and on fur- 
ther re moval of unbound galaxies id entified in projected phase- 
space dden Hartog & KatgertI 1 19961) . a procedure validated via 
the c omparison wit h cluster s extracted from numerical simula- 
tions (IBiviano et al . 2006; W oitak et al.ll2007l) . On the remaining 
cluster members we determine cluster los velocity dispersions 
cT/oj using the robust biweight estimator (iBeers et al. 199(]j). 
Finally, we determine M20 cluster masses adopting the scal- 
ing 0-10S-M20Q relation of iMauduit & Marnoiil d2007 ', app. A), 
which is in go od agreement wi th the phenomenological relation 
derived by Biviano et al. (2006) for a sample of cluster-sized ha- 
los extracted from cosmological simulations. 

In summary, our data-set consists of 19 distant clusters from 
z = 0.393 to 0.794 with a mean (median) redshift z = 0.56(0.54), 
and 59 nearby clusters from z - 0.035 to 0.098 with a mean 
(median) redshift z = 0.07(0.06). The distant clusters span the 
M200 mass-range 0.7 - 13.6 x IO'^^Mq, with a mean (me- 
dian) mass of 2.8(4.4) x 10^* Mq. The nearby clusters span the 
M200 mass-range 0.4 - 20.5 x lO''* Mq, with a mean (median) 
mass of 5.9(5.7) x 10'"*Mq. Assu ming a cluster mas s accretion 
rate of ~ 0.08M2oo(z = 0)/Gyr dAdami et alj|2005h . the aver- 
age mass of our high-z cluster sample is expected to increase to 
about the average mass of our low-z cluster sample in the time 
that separates the two cosmic epochs (~ 5 Gyr). A similar, albeit 
some what larger, evolution in mass is also predicted theoreti- 
cally dLapi & CavaUerell2009l) . In this sense, we are comparing 
similar objects observed at two different cosmic epochs. 

The low-z and high-z cluster samples are presented in Table[T| 
and|2l respectively. In col.(l) we Ust the cluster n ame (follow- 
ing the short-name convention of iPoggianti et alj|2006. for the 
high-z clusters), in col. (2) its mean redshift, in col.(3) the cluster 
M200 in 10^* Mq units, in col. (4) the number of galaxies without 



' The virial radius r2oo is the radius within which the enclosed aver- 
age mass density of a cluster is 200 times the critical density. 
The virial mass M200 is the mass enclosed within a sphere of radius 

'■200- 

The circular velocity is defined from the previous two quantities as 

V2oo=(G Mzoo/f-^oo)"^- 
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emission lines in their spectra (nELGs, hereafter) in the radial Table 1. The low-z cluster sample, 
range 0.05 < /?/r2oo ^ 1, in col. (5) the number of galaxies with 
emission lines (ELGs, hereafter) in the same radial range, and in 
col. (6) the radial distance from the cluster center of the most dis- 
tant galaxy in the sample, in units of r2oo. The ELG classification 
for the low-z cluster sample is described in Katgert et al. ( 199^. 
For the high-z cluster sample we classify ELGs the EDisCS 
galaxies with an [Oil] equivale nt width > 3 A or with any 
other line in emission (Poggiant i et al.ll2006l) and the MORPHS 
galaxies with a spectra l type different from 'k', 'k+a', and 'a+k' 
jPoggianti et~anil999l) . 

3. The construction of the stacked cluster samples 

In order to be able to analyse the cluster mass and velocity 
anisotropy profiles, the available spectroscopic data of individual 
clusters from our data-sets are not sufficient. We need to stack 
all the clusters from each of the two samples together. Stacked 
cluster samples have been used succes sfully in several analy- 
ses of the properties of clusters (e.g. lMoss~& Die ken s| 1 19771; 
Biviano et al.l 119921: IC arlberg et al? 'I997at lv~n der Marel et alj 
2000; Katgert etalJE004; Rines et al. 200%. The validity of this 
approach is supported by the results of cosmological numerical 
simulations that predict cosmological halo s to be characterize d 
by the same, universal mass density profile jNavarro et al.119971) . 
Even if the mass density profi les of cosmological halos do de- 
pend on their mass (see, e.g., iNavarro et al. 1 1 19971: iDolag et alj 
12004 ). this dependence is very mild, and the profiles are very 
similar for halos with masses within about two decades around 
the average cluster-like halo mass. In order to ensure homology 
of our distant cluster data-set we rejected three very low-mass 
clusters from the initial sample (see Sect.|2]i. 

The observables on which the analysis is based (see 
Sect, m are the galaxy projected clustercentric distances, R, 
and the galaxy los velo cities in their cluster rest frames 
jHarrison & Noonanlll979l) . \rf = (v - v7)/(l + z^, where v 
are the observed galaxy velocities, and z^, v^ the average cluster 
redshift and los velocity, respectively. For the cluster centers we 
use the position of the X-ray surface brightness peak, or, if this 
is unavailable, the position of the cluster brightest cluster galaxy 
(BCG hereafter). Since different determinations of a clu ster cen - 
ter typically diff'er by less than 100 kpc ( Adami et"anil998bh . 
precise cent ering is not very imp ortant for our dynamical analy- 
sis (see also lBiviano et al.l l2006^. 

In order to stack clusters together it is necessary to adopt ap- 
propriate cluster-dependent scalings for the R and Vr/ quantities. 
Stacking the clusters in physical units would have the unwanted 
consequence of mixing up the virialized regions of some clus- 
ters with the unvirialized, external regions of others. Radii and 
velocities are therefore scaled by the clusters virial radii, r2oo, 
and circular velocities, V20(), so that the dynamical analysis on 
the stacked cluster is done in the normalized units R,, = R/r2oo 
and v„ = Vr//v2oo- 

In our analysis we consider only the virialized clus- 
ter regions, i.e. galaxies with R„ < I. The Jeans method 
dBinnev & Tremaind 19871) that we adopt here to determine clus- 
ter mass and anisotropy profiles is in fact not valid outside the 
region where dynamical equilibrium is likely to hold. Moreover, 
we also exclude the very central cluster regions (R„ < 0.05) in 
order to account for the positional uncertainties of cluster cen- 
ters, and in order not to include the centrally located BCGs 
in our sample. BCGs are probably built up via merger pro- 
cesses even quite recently (see, e.g., De Lucia & Blaizoj|2007t 
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iRamella et al.ll2007t iRines et al.ll2007h . Including these galaxies 



would therefore probably invalidate t he Jeans analysis in which 
dissipation processes are not included dMenci & Fusco-Femiand 
,1996) . 

In the selected radial range 0.05 < R„ < I there are 556 
galaxies in the distant stacked cluster, and 2566 in the nearby 
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Table 2. The high-z cluster sample. 
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one. We further split these samples in two by considering nELGs 
and ELGs separately. There are 293 nELGs and 263 ELGs in 
the distant stacked cluster, and 2226 nELGs and 340 ELGs in 
the nearby one. The much larger fraction of ELGs in the dis- 
tant stacked cluster compared to the nearby one is a known fea- 
ture of the evolutioii of the cluster galaxy popu lation fractions 
(iDressler et al.lll999l:IPoggianti et all 1999112006.) . 

Neither the low-z nor the high-z cluster samples are spectro- 
scopically complete to a given magnitude. Incompleteness is not 
a problem in the dynamical analysis as far as it is the same at 
all radii. In fact, the normalization of the galaxy number den- 
sity profile cancels out in the Jeans procedure (see eqs.lO and 
(IU in Sect. HJl and the shape of the galaxy number-density and 
velocity-dispersion profiles are known to be independent from 
the galaxy luminositi es, at least for absolute magnitudes Mr > 
-22.8 dBiviano et al.ll2002l) . Galaxies brighter than Mr = -22.8 
are mostly BCGs, and have been mostly excluded from our sam- 
ples by removing the very central cluster regions, R„ < 0.05. 
While there are several indications that dwarf galaxies in clus- 
ters hav e a differe nt phase-s pace dist ribution from bright galax- 
ies (e.g..lLobo et a l. 1997; Kambas et al. 2000; Odell et al. 2002; 
iPopesso et al.ll2006 ), this is irrelevant here since dwarf galaxies 
are not present in our samples. 

If incompleteness does depend on radius, a correction must 
be applied. In fact, a radially-dependent incompleteness would 
change not only the normalization of the galaxy density profile 
but also its shape. For the ENACS, from which our local sam- 
ple is extracted, it has been shown t hat the spectrosco pic incom- 
pleteness is not radial dependent (iKatgert et al.lll998l) . On the 
other hand, for the high-z cluster sample spectroscopic incom- 
pleteness is a function of radius, althoiigh no t a strong one, and 
must be connected for ("Poggian ti et al.ll2b06h . The incomplete- 
ness correction for each cluster in our high-z sample is done by 
applying the geometri cal weights with the method described in 
iPoggianti et al.l (l2006i Appendix A). 

Another kind of radial incompleteness results from the stack- 
ing procedure when the clusters that enter the stacked sample do 
not cover the same radial range, i.e. they have not all been sam- 
pled out to the same limiting aperture. This happens for both our 
local and distant samples (the apertures are listed in the last col- 
umn of Tables [T] and IU. A radical way of addressing this prob- 
lem would be to stack clusters only out to the minimum among 



Fig. 1. The projected number density profiles, N(R„), of the 
nELGs and ELGs (upper and lower panel, respectively) for the 
low- and high-z stacked clusters (left and right panel, respec- 
tively). Solid lines represent best-fit models to the data, i.e. pro- 
jected NEW models for nELGs and high-z ELGs, and the core 
model for low-z ELGs (see Sect. |5.1| and |5.2| i. 



all clusters apertures. In order to make the most efficient use of 
our sample, and explore the clusters dynamics out to r2oo, we 
prefer to follow another approach. We apply a correction for the 
fact that not all clusters contribute at all radii. This coiTection is 
based on the relative individual cluster contributions to the to- 
tal number of galaxies in the stacked sample, estimated at the 
radii where the clusters are still sampled. E.g. if at a given ra- 
dius R„^i„, there are m clusters that have not been sampled out to 
that radius, their contributions at radii R„ > R„ ,„ is accounted 
for by applying a correction factor 1/(1 - 2™ 1 fd to the total 
galaxy counts in the stacked cluster, where ft is the fraction of 
galaxies contributed to the total sample by cluster / in the radial 
range where the cluster is sampled. The m ethod is similar to the 
one described in lMerrifi eld & K eli3(ll989l) an d has been applied 
to the ENACS sample in Bivian o'etal.l j2002l) and iKatgert et all 
(^2004). 

In Fig.[T]and|2]we display the projected number density pro- 
files, N{Rn), and the los velocity dispersion profiles, criosiRn), re- 
spectively, of nELGs and ELGs for our local and distant stacked 
clusters. 



4. The dynamical analysis: method 

The method we adopt for the dynamical analysis of our two 
stacked cluster s is based on the standard spherically-symmetric 
Jeans analysis dBinnev & Tremain3ll987h . 

The assumption of spherical symmetry appears justified by 
the fact that our samples combine ma ny clusters together, irre- 
spective of their or ientation (see also van der Marel et aLll2000l; 
IKatgert et al.ll2004 . The resulting stacked clusters are spheri- 
cally symmetric by construction, except if the cluster selection 
process favors a preferential orientation. E.g., it has been argued 
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Fig. 2. The los velocity dispersion profiles crios(Rn), of the 
nELGs and ELGs (upper and lower panel, respectively) for the 
low- and high-z stacked cluster (left and right panel, respec- 
tively). Solid lines represent best-fit models to the data, i.e. 
c = 4.0 NFW mass model, a = 3.6 (respectively a = 1.2) 
OM velocity-anisotropy model for low-z nELGs (respectively, 
ELGs), and c = 3.2 NFW mass model, a = 0.01 ML velocity- 
anisotropy model for high-z nELGs and ELGs (see Sect. 15. H and 

ED- 

that clusters selected because of the presence of gravitational 
arcs are more likely to have their major axes aligned along the 
los (e.g. lAlleal998) . Neither ENACS nor EDisCS clusters were 
selected because of the presence of gravitational arcs, and they 
are both likely to be representative of the overall cluster pop- 
ulations in their mass and redshift range s. ENACS clusters are 
drawn from the catalog of rich clusters in lAbell et al.l(fl989 ). The 
mass function derived using clusters from this ca talog is simi- 
lar to that obtained using X-ray sele cted clusters dMazure et al.l 
11996 ^. hence rich clusters from the Abel l et al. I (119891) catalog 
appear to be representative of all rich clusters in the nearby uni- 
verse. EDisCS clusters were optically selected from the high- 
est surface brig htness candidates in t he Las Campanas Distant 
Cluster Survey (iGonzalez et al.ll200TI) . From a comparis on with 
other published optical/X-ray cluster catalogs, Gonzal ez et al.l 
(1200 Ih have shown that their detection method is able to recover 
> 90% of the cluster population. Further support against a pos- 
sible orientation bias comes from the comparison of different 
EDisCS cluster mass-estimates, obtained using the distribution 
of cluster galaxies, gravitat ionally lensed images (Clo we et al.l 
120061: iMilvang- Jensen et al.'2008 ) and the X -ray emission from 
the infra-cluster plasma (Johnso n et al.ll20d6l) . 

In the spherically-symmetric Jeans analysis, our observables 
are the galaxy number density profile N{R„) and the los velocity 
dispersion profile cr;„j(/?„). N(Rn) is uniquely related to the 3- 
dimensional (3-d) galaxy number density profile v(r„) via the 
Abel inversion equation. 



1 r 

v(r„) = — 



dN dR„ 



dRn V^F^ 



(1) 



where r„ = rlrioa is the 3-d clustercentric radius in normalized 
units. The other observable, a-i„s{R„) is related to the cluster mass 
profile, M(r„), and the cluster velocity anisotropy profile. 



/^(r„) = 1 



<v?> 



(2) 



where <v^>, <v^> are the mean squared tangential and radial 
velocity components, which reduce to cr^ and cr^ respectively, in 
the absence of bulk motions and net rotation. Given M(r„) and 
'(r„) the observable criosiRn) follows through (.van der Marel 
1991 



I 



, f"" vM 

o-,.{rn) = G/v(r„) J — exp 

and dBinnev & Mamonll 19821) 



f3dx 

X 



d^. 



(rfJRn) = 2/N(R, 



«) r 



Rt \ vcT].r„dr„ 



(3) 



(4) 



where G is the gravitational constant. Note that in practice the 
upper limit of the integrals in the above equations is set to a 
finite radius (typically 20 in the integration units), large enough 
as to ensure that the result of the integration does not change 
significantly by pushing the limit to larger values. 

It is therefore possible to adopt parameterized model repre- 
sentations of M{r„) and y6(r„) and determine the best-fit param- 
eters by comparing the observed crigs{R„) profile with the pre- 
dicted one, using the statistics and the uncertainties on the 
observed profile. 

From the eqs. above it is however clear that different com- 
binations of the mass and anisotropy profiles can produce 
the same los velocity dispersion profile, the so-called "mass- 
anisotropy" degeneracy. Different methods exist to solve this 
degeneracy (see, e.g., Merritt 1987; van der Marel et al. 200^ 
iLokas & Mamonll2003l: IWu & Tremainell2006i) . These methods 
are effective for data-sets of ~ 1000 tracers of the gravita- 
tional potential. Given the smaller size of our distant clusters 
data-set we here adopt another method, recently suggested by 
Battaglia et al. (2008) for application to the case of dwarf galax- 
ies. 

The method of lBattaglia et al.l (l2008 l) consists in considering 
not one, but two different tracers of the cluster gravitational po- 
tential, so that there are two observables (the los velocity disper- 
sion profiles of the two tracers) to solve for the two unknowns, 
M{r„) and /3{r„). Since M(r„) must be the same for both trac- 
ers, but y6(r„) can in principle be different, the degeneracy is only 
partially broken, however the constraints on the dynamics of the 
system are significantly stronger than with a single tracer Of 
course, this method works only if the two tracers have different 
projected phase-space distributions. This is the case of cluster 
nELGs and ELGs, which we therefore adopt as our two popula- 
tions of tracers. 

In practice, we adopt parameterized models for M{r„) and 
/3(r„), solve eqs.([TJ, (O, and (|4]l separately for nELGs and ELGs 
in each of the two stacked clusters, and jointly compare the so- 
lutions to the observed criosiRn) of the two galaxy populations. 
Say ctoIjsj and cr,„odei,i the observed and predicted crios of a given 
population of tracers (nELG or ELG), in the /-th of m radial bins, 
and say 5, the corresponding l-cr uncertainty on cTohsj- We find 
the best-fit parameter of the model representing M(r„) and its 
uncertainties by minimizing 



2 _ 2 ,2 
X - XnEW XeW 



(5) 
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Xtracer 



(6) 



Since M{r„) must be unique for the two populations, we then 
adopt the joint best-fit solution for M(r„) and solve again eqs.lO 
and (HI separately for nELGs and ELGs to determine the best-fit 
parameters (and uncertainties) of their j6-profiles through the 
minimization, eq.©, i.e. we marginalize over M{r„) to constrain 
the two y6(r„) solutions. 

The choice of the models cannot be too generic nor too re- 
strictive. If it is too restrictive, we might find it difficult to find 
accurate fits to the data. On the other hand, with small data-sets 
(such as that of our high-z sample) it would be difficult to obtain 
strong constraints on models that are too generic or are charac- 
terized by too many parameters. We let ourselves be guided in 
our choices by theoretical expectations. 

The distribution of mass within cosmological halos, such as 
galaxy clusters, is a robust prediction of the Cold Dark Matter 
(CDM) cosmological model. Numerical simulations have shown 
that all cosmological halos are characterized by the same, uni- 
versal, mass density profile (.Navarro et al.iil997i) . the so-called 
'NFW profile, 

Pnfw °^ (cr„y\l + cr„y^, (7) 

characterized by the concentration parameter c, a central cusp, 
and an asymptotic slope of -3 at large radii. While other 
analytical forms have subsequeiitly been proposed (see, e.g . 
iMoore et alj 119991: iHavashi et al] 120041: iDiemand et al] l2005i) . 
the NFW profile provides an acceptable fit to observed mass pro- 
files of ga laxy clusters both at low and intermediate redshifts 



(see, e.g. Ivan der Marel et alj l2000t iBiviano & Girardil l2003t 



iKatgert et al.ll2004l:lz^pacosta et alj2006l) . It is therefore rather 
straightforward to choose the NFW profile as our reference 
model for M(r). 

We adopt the NFW profile also as a model for v( r„) (actually, 
we fi t N(R„) with the projec ted NFW profile, see iBartelmantj 
[T99llLokas & Mamonll200"l . but we do not make the assump- 
tion that v(r„) and M(r„) are characterized by the same NFW 
profile. I.e. we do not work in the so-called light-traces-mass hy- 
pothesis. Whenever the NFW model does not provide an accept- 
able fit to v(r„), we consider a different model, with one addi- 
tional parameter, to allow for a bett er fit. The model we adopt in 
this c ase is the so-called j6-model jCavaliere & Fusco-Femianol 
fT978h . 



N^il+iRnlR.fr" 



(8) 



characterized by two parameters, the core-radius, Rc, and the 
slope, a. In order to avoid terminology confusion with the 
velocity-anisotropy profile y6(r„), in the following we refer to the 
j8-model as the 'core' model. 

The model choice for the velocity-anisotropy profile j6(r„) is 
less straightforward. At variance with the case of the halo mass 
profile, there is no claimed 'universal' velocity anisotropy pro- 
file. We consider the two follo wing models. One is the Mamon- 
Lokas ('ML' hereafter) model dMamon & Lokasll2005l) 



jS = 0.5 r„/(r„ + a). 



(9) 



which has been shown to provide a good fit to the velocity 
anisotropy profiles of simulated cosmological h alos. The other 
is the Osipko v-MeiTitt ('OM' hereafter) model jOsipkovll 19791: 
lMerrittlll985h 



which provides a good fit to the observed velocity 
anisotropy profile of late -type galaxies in the ENACS sample 
(iBiviano & KatgertI l2004h . Both the ML and the OM models 
depend on just one free parameter, the anisotropy radius, a. 



5. Results 

5.1. The nearby cluster sample 

The NFW profile (in projection) provides an acceptable fit to the 
nELG N(R„), with c = 2.4. This best-fit profile is shown in the 
upper-left panel of Fig.[T] On the other hand, the ELG N{R) can- 
not be fitted by a (projected) NFW profile, because the ELGs 
avoid the central cluster region. The wider spatial distribution 
of the ELGs as compared to that of the nELGs clearly reflects 
the well known MDR (see Sect. [T]). In order to fit the rather flat 
ELG N(R„), we adopt the core model, and we find an accept- 
able fit with Re - 1.28 and a - -3.2 (see the bottom-left panel 
of Fig. [T]i. The N(R„) best-fitting models are Abel-inverted to 
provide the 3-d number density profiles v(r„) that we use in the 
Jeans analysis. 

The observed cr/oj-profiles of the nELGs and ELGs are 
shown in the upper-left and, respectively, lower-left panels of 
Fig.|2] The strikingly different behavior of the two cr/o, -profiles 
reflects the well known morphological segregation in velocity 
space (see Sect.[TJ. We then determine the best-fit concentration 
parameter of a NFW M{r„) model by a joint;^^^ comparison of 
the observed nELG and ELG cr/o^-profiles with those predicted 
by the Jeans analysis for the given mass model, and for either 
ML or OM j6(r„) models. The best-fit solution is obtained using 
the OM model for the velocity anisotropy profiles of nELGs and 
ELGs. The best-fit value of the NFW concentration parameter 



is c = 4.0 



(+0.5)+2.3 
(-0.1)-1.3 



(90% confidence levels, c.l. in the following. 



(10) 



68% c.l. in brackets). The;^^^ vs. c solution is displayed in the 
top panel of Fig. [3] 

Marginalizing over this Mi/n) best-fit solution we obtain the 
best-fit y6(r„) OM-model parameters a - 3.6^_^^^ (in units of r2oo) 
and a - \ .2^q^ for the nELG and ELG populations, respectively. 
We only provide here 90% c.l., since we find no acceptable so- 
lutions at the 68% c.l. for the ELGs. The best-fit model cr/o^- 
profiles are shown as solid lines overlaid on the observed, binned 
ciosiRn) in Fig. 12] (left-hand panels). Clearly, the fit to the ELG 
(TiosiRn) is not excellent, but still acceptable to within 90% c.l. 

The solutions for the velocity-anisotropy profiles are dis- 
played in Fig. |4] for the nELGs (upper-left panel) and the ELGs 
(lower-left panel). Notice that the quantity displayed in Fig. |4] 
is (Tr/cr, = (1 -yS)"'^^. The nELG velocity-anisotropy profile is 
consistent with fully isotropic orbits of this population of galax- 
ies withi n the virial region, i n agreement with the results ob- 
tained bv iKatpert et all ('2004). On the other hand, the ELG or- 
bits are approximately isotropic only out to r„ ^ 0.6, and then 
become increasingly radial. The OM-model solution is very sim- 
ilar to the non-parametric velocity-anisotropy profile derived by 
IBiviano & Katgert (2004). It is also quite similar to the velocity- 
anisotropy profi les of DM particles in cluster-size simula ted ha- 
los at z » (e.g. lTormen et al ]|l997HDiaferio et al.ll200ll) . 

For the sake of comparison with the distant cluster sample 
(see Sect. I5.2l i it is useful to also consider the solution obtained 
for the ML y6(r„) model. With such a model, we obtain an ac- 
ceptable solution of the Jeans analysis for a NFW mass profile 
with the same concentration value obtained using the OM I3{r„) 
model (c = 4), although with a larger value. Marginalizing 
over the c - A parameter, we then obtain the best-fit yS(r„) ML 
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Fig. 3. The values vs. the c parameter of the NFW mass pro- 
file model for our low-z (top panel) and high-z (bottom panel) 
stacked clusters. The two dashed lines indicate the 68% and 90% 
C.I., and the vertical lines indicate the best-fit c values. 
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Fig. 4. Best-fit velocity-anisotropy profile (Tr/cT, = (1 -y6)"'^^ for 
nELGs (upper left panel) and ELGs (lower left panel) in low- 
z clusters and for nELGs (upper right panel) and ELGs (lower 
right panel) in high-z clusters. The best-fit velocity-anisotropy 
model is OM for the low-z sample and ML for the high-z sample. 
Dashed lines indicate the 90% c.l. For the high-z sample, the 
best-fit solutions are at the lower-limit of the interval considered 
in the;if^ minimization analysis, hence only the lower c.l. to /3{r„) 
Sire shown. 



model parameters a = 2.9lf^ and a = 1.5^^^ (90% c.l.) for the 
nELG and ELG populations, respectively. Similarly to what is 
obtained using the OM y6(r„) model, also in this case the best-fit 
anisotropy radius of the nELG population is about twice as large 
as that of the ELG population, indicating that the ELG orbits in 
clusters are more radially elongated than the nELG orbits. 



5.2. The distant cluster sample 

Both the nELG and the ELG N{R„) are well fitted by projected 
NFW profiles with c = 7.5 and c - 2.7, respectively. These 
N(R„) best-fitting models are Abel-inverted to provide the 3-d 
number density profiles v(r„) that we use in the Jeans analysis. 
The binned profiles and their best-fit models are shown in the 
right-hand panels of Fig. [1] Similarly to what is found in the 
nearby cluster sample, ELGs have a wider spatial distribution 
than nELGs, confirming previous results about the existence of 
a MDR also at z ~ 0.6 (see Sect.[TJ. However, the ELGs do not 
avoid the central cluster regions as in low-z clusters. 

The evidence for segregation in velocity space is not so 
strong. In fact the cr/„5-pro files of nELGs and ELGs are not very 
different (see the right-hand panels of Fig. |2]l, in contrast with 
the situation seen in low-z clusters. We jointly compare these 
cr/o5-profiles with those predicted by NFW M(r„)-models, and 
either ML or OM /3{r„) models, via the method described 
in Sect. |4] The best-fit solution is obtained for the ML /3(r„) 
model. The best-fit value of the NFW concentration parameter 
is c = 3.2|^J o*:^2o (90% C.L, 68% c.l. in brackets). The vs. c 
solution is displayed in Fig. |3] bottom panel. As expected from 
the smaller size of the high-z sample, the solution is less well 
constrained than for the low-z sample. 

We then adopt the c = 3.2 NFW best- fit solution as the refer- 
ence M{r„) for the stacked high-z cluster, and look for the best- 
fit y6(r„) ML-model solutions, separately for nELGs and ELGs. 
We find that the best-fit anisotropy-radius parameter is identical 
for the two galaxy classes, a - 0.01, at the lower limit of the a- 
range considered in the;^^^ minimization analysis. Unfortunately, 
the solutions are poorly constrained. The 90% upper limit for the 
nELGs is a < 0.9, that for the ELGs is at the upper limit of the a- 
range considered, a < 10.0. The best-fit model cr/g, -profiles are 
shown as solid lines overlaid on the observed, binned a-ioAfn) in 
Fig. 12] (right-hand panels). 

The solutions for the velocity-anisotropy profiles, 0-^/0-,, are 
displayed in Fig. |4] (right-hand panels) and are identical for the 
nELGs and the ELGs. The upper limits on the a-parameters 
translate in lower-envelopes to the velocity-anisotropy profiles 
(dashed curves in the right-hand panels of Fig.|4]i. Taken at face 
value, both the nELGs and the ELGs have radially anisotropic 
orbits. Isotropic orbits are excluded for high-z nELGs, at least 
outside the center, but cannot be excluded for high-z ELGs be- 
cause of the large error bars. 

The poor constraints on the velocity anisotropy of high- 
z ELGs do not allow us to draw any conclusion about their 
orbital evolution with z. Taken at face value the results sug- 
gest that ELGs follow radially anisotropic orbits both in high-z 
and in low-z clusters. On the other hand, the orbits of nELGs 
evolve with z, from almost isotropic in low-z clusters to radially 
anisotropic in high-z clusters (compare the left-hand and right- 
hand panels of Fig. |4]i. 
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Fig. 5. The distributions of predicted c for the low-z (upper 
panel) and high-z (l ower pane l) clust er samples, as obtained 
from the re l ation of iGao et al] (l2008l) (solid histograms) and 
iDuffv et alj (|2008|) (dashed histograms) applied to our cluster 
mass and redshift estimates. The dots indicate the best-fit c val- 
ues from the Jeans analysis, and the solid and dashed lines above 
the histograms the 68%, respectively 90%, confidence intervals 
on the best-fit solutions. 



6. Summary and discussion 

We have obtained acceptable equilibrium solutions to the nELG 
and ELG projected phase-space distributions with NFW M{r„) 
models, and with OM and, respectively, ML yS(r„) models, for 
the low-z and, respectively, the high-z sample. The solutions do 
not always provide excellent fits to the crios{Rn) of the nELGs 
and ELGs, but they are anyway always acceptable at the 90% 
c.l. Hence, the rather limited set of models we have considered 
seems to be adequate for our data-sets. Larger data-sets would 
be required if a wider range of models and/or parameters is to be 
considered. 

The best-fit NFW c values found with the Jeans analysis are 
in agreement with those expected in a ACDM universe for cos- 
mological halos with the masses and redshifts of the clusters in 
our samples. In Figure |5] we display the distribution of the pre- 
dicted c values for the clusters in the low-z and high-z samples. 
The predicted c values are obtained from our cluster mass and 
redsh ift esti mates, usin g the c = c(M,z) relations of iGao et al.l 
(I2OO81) and iDuff'v etal.l (|2008|) . In the same figure we also dis- 
play the best-fit c value obtained with the Jeans analysis, and the 
68 and 90% c.l. The Jeans solutions are fully consistent with the 
predicted distributions of c values, both at low- and at high-z. 

iDuffv et aTl (l2008h claimed that a discrepancy exists be- 
tween the predicted c values for groups and clusters of galaxies 
and those determined observationally using X-ray observations. 
Since this is not apparent here, the claimed dis crepancy may oc- 
cur mostly at the group scale (see Fig. 8 in iBivianol 12008) or 
originate from a systematic bias in the the X-ray-based cluster 
mass estimates dRasia et al.ll2006l) . 



Our results are therefore consistent with the theoretical c - 
c(M, z) relation for cosmological cluster-sized halos in a ACDM 
universe. 

As far as the orbits of cluster galaxies are concerned, our 
analysis suggests that these orbits become more isotropic with 
time. Low-z cluster nELGs have nearly isotropic orbits, low-z 
cluster ELGs have radially anisotropic orbits outside the central 
cluster regions. We can exclude that ELGs follow isotropic orbits 
at the 90% c.l. outside r ^ 0.6 rioo. 

In high-z clusters nELGs and ELGs appear to be character- 
ized by similar orbital anisotropics, radial and almost constant 
with radius. For nELGs we can exclude isotropic orbits outside 
r ^ 0.1 r2oo at the 90% c.l. Since low-z cluster nELGs are char- 
acterized by nearly isotropic orbits, the orbits of nELGs must 
evolve from radial to isotropic with time. No significant evo- 
lution is found for the orbits of ELGs, which remain radially 
anisotropic both at high- and low-z. 

A process capable of isotropizing galaxy orbits with time 
is the secu lar growth of cluster mass via hierarchical accretion 
(iGill et al.l 120041) . According to theoretical models, the growth 
of halo masses occurs in two phases, an initial, fast accre- 
tion phase, followed by a slower, smoother accretion phase 
(iLapi & Cavalierell200^ . During the fast accretion phase clus- 
ters undergo major mergers that can induce rapid changes 
in the cluster gravitational potentia l jManrique et al.l 120031; 
iPeirani et al.ll20^ IValluri et al.ll2007h . causing energy and an- 
gular momentum mixing in the galaxy distributions, and thereby 
isotropization of galaxy orbits (Henon 1964; Lvnden-Bell 1967f 
Kandrup & Siopis 2003; Merritt 2005; Lapi & Cavaliere 2009). 

Cluster-sized halos of 10''*~'^Mo mass undergo the transi- 
tion fro m the fast to the slow a ccretion phase at z ~ 0.4 (see 
Fig.6 in iLapi & Cavaiiergi2009l) . Hence, clusters in our high-z 
sample are observed before the end of their fast accretion phase. 
During this phase, the orbits of cluster galaxies can still evolve, 
approaching isotropy. At z < 0.4 one therefore expects to see 
isotropic orbits for those galaxies that were already part of the 
clusters before the end of the fast accretion phase. In our low-z 
cluster sample we therefore expect nELGs to be characterized by 
isotropic orbits, as observed. On the other hand, ELGs must be 
newcomers in low-z clusters, hence memory of their recent clus- 
ter infall is still conserved in their (radially elongated) or bits. 

The fact that a significant fraction (23 + 3%, Poggiant i et al.l 
2009) of the high-z nELGs have ak + a 01 a + k spectral classifi- 
cation is supporting evidence for the fact that the clusters of our 
high-z sample are still in the fast accretion phase. The spectral 
characteristics of these galaxies are in fact indicative that their 
star formation have stoppe d in the ~ 1 .5 Gyr prior to observa- 
tion jPoggianti et al.|[T999l) . If the cessation of the star-formation 
activity in these galaxies is related to their first encounter with 
the cluster environment, their accretion has been relatively re- 
cent. Hence, a large part of the cluster galaxies (and presumably, 
of the cluster mass) has been assembled in the last ~ 1 .5 Gyr. 

The detected orbital evolution occurs over the same period 
of cosmic time (the last ~ 5 Gyr) when the cluster galaxy 
population undergoes a major change in its mo rphological mix 
(De sai et al.l2007l) and star formation properties dPoggianti et al.l 
2006^ We might then be witnessing two effects of the same un- 
derlying physical phenomenon. 

The transformation of ELGs into nELGs may be at least 
partly responsible for another evolutionary trend we observe in 
our clusters, that of the nELG number density profile. This pro- 
file becomes less concentrated with time, an evolution in the op- 
posite sense to that observed for the mass density profile. The 
(projected) NFW models fitted to the number density profile of 
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low- and high-z nELGs have best-fit concentrations c - 2.4 ^Jj j 
and c - 7.5^y g (90% c.l.; see also Figure[TJ. On the other hand, 
no significant evolution is found for the number density pro- 
file of nELGs and ELGs together. This profile is dominated by 
nELGs at low-z, but not at high-z. If high-z ELGs transform into 
nELGs with time, the nELG-i-ELG number density profile would 
not change, but the nELG number density profile would flatten, 
since ELGs are less spatially concentrated than nELGs. Also the 
ELG number density profile flattens with time. This might be re- 
lated to the cluster environment growing more hostile with time, 
and making more difficult for infalling field galaxies to conserve 
their gas as they approach and cross the cluster centers. 

The results we have obtained in the present study are based 
on the still rather limited amount of available data for high-z 
cluster galaxies. Moreover, our high-z and low-z cluster samples 
span quite a substantial range in masses. It will then be impor- 
tant to tighten the current constraints on the orbital evolution of 
cluster galaxies using future, larger spectroscopic data-sets for 
high-z clusters, and also to re-assess such an evolution as a func- 
tion of cluster mass. From these future analyses we will obtain 
a more thorough understanding of the hierarchical assembly his- 
tory and evolution of galaxy clusters. 
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